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AN EXPERIMENTAL AND NUMERICAL STUDY OF 
WAVE MOTION AND UPSTREAM INFLUENCE 
IN A STRATIFIED FLUID 

David A. Hurdis 
Goddard Space Flight Center 


INTRODUCTION 

Wave Propagation in Nonhomogeneous Fluids 

The earliest significant study of waves in a nonhomogeneous fluid was that of Stokes (1847) 
who developed the theory of the propagation of infinitesimal waves in a two-layer fluid 
system. The excitation of waves in two fluid layers is discussed also by Lamb (1932, 

Section 231). In essence, it has been shown that for a two-layer fluid, two possible wave 
systems can exist for any given period of oscillation. One propagates at the free surface 
of the upper layer, while the other propagates at the interface of the two-layer fluid. 
Similarly, it can be shown for a three-layer fluid that three wave systems can exist, one for 
each surface of density discontinuity. In this context, a fluid with a continuous density 
stratification can be considered to be the limit of an n-layer fluid system as n -»■ internal 
waves can be propagated at all levels in such a fluid. A semantic distinction is usually made 
between waves propagating at the interfaces of a multilayer system and those propagating 
through a continuously stratified fluid. The former are called interfacial waves, while the 
latter are referred to as internal waves. 

Internal wave propagation in a continuously stratified fluid is characterized by the local 
Brunt-Vaisala frequency N = (-g/p dp/dz) 1/2 . (For stably stratified fluids, dp/dz < 0.) No 
waves can propagate with frequencies higher than the local value of N, and the direction of 
propagation of a wave of frequency a is cos' 1 (o/N) from the vertical. 

The earliest important experimental work on wave propagation in a nonhomogeneous fluid 
was published by Ekman (1904), who sought to explain the abnormally high resistance 
experienced by ships in the mouths of some Norwegian fjords where there is a layer of 
fresh water over salt water. Ekman demonstrated that the phenomenon is due to increased 
wave drag resulting from the ship’s exciting of waves not only at the free surface but also 
at the interface of the two layers. 

Upstream Influence 

Density stratification affects the flow of an inviscid, incompressible fluid in two ways. 

First, a change in density causes a change in the inertia per unit volume of the fluid. 

Second, in a gravitational field, a density change causes a change in. the body force per 



unit volume of the fluid. Since, in general, the body force is not the only force acting on 
the fluid, these two effects will not nullify each other, and the flow pattern of a nonhomo- 
geneous fluid may differ considerably from that of a homogeneous fluid even though the 
boundary conditions are the same. 

One flow phenomenon exhibited by stratified fluids is upstream influence. This phenome- 
non occurs during the flow of a stratified fluid past an obstacle, and consists of the propa- 
gation of a disturbance upstream from the obstacle. It is especially strong when the obsta- 
cle is so wide, in comparison with its height, that the flow past the obstacle is essentially 
two-dimensional. Upstream influence is also called forward influence or upstream wake in 
the literature. The limiting case and strongest manifestation of upstream influence occurs 
in the slow flow of a continuously stratified fluid past an obstacle, and is called blocking. 
Blocking is the stagnation with respect to an obstacle of the fluid within a horizontal slab 
containing the obstacle and having the same height as the obstacle. Gosed streamline 
recirculation may occur within the slug of blocked fluid being pushed ahead of the obstacle 
(Janowitz 1971 ; Browand and Winant 1972). For very low flow speeds, blocking occurs 
both upstream and downstream of the obstacle. At higher speeds, the downstream blocking 
is replaced by a wake and lee waves, and as the speed increases still more, the upstream 
blocking becomes incomplete and is therefore more properly referred to as upstream influ- 
ence. 

Qualitative explanations for upstream influence may be given from two viewpoints, the 
buoyancy viewpoint and the internal wave viewpoint. 

Buoyancy Viewpoint 

As a two-dimensional body moves through a stably stratified fluid, it is necessary for the 
fluid to be displaced vertically to get by the body. Because of the stable stratification, 
however, the fluid is reluctant to be displaced vertically, and tends to pile up in front of 
the body, thus accounting for the upstream influence. 

Internal Wave Viewpoint 

The upstream influence may also be regarded as the propagation of an internal wave of 
zero frequency. From the above discussion of internal waves, it is clear that such a wave 
can propagate from the disturbance source only in a horizontal direction, and that as the 
wave does so, it causes the presence of the disturbance to be felt far upstream. 

Investigations of Upstream Influence 

Experimental 

The fundamental experiment which established the existence of the phenomenon of 
upstream blocking in stratified flows was conducted by Long (1955). He studied the flow 
of a linearly stratified salt solution over a two-dimensional barrier which was towed along 
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the bottom of a channel. He observed that for values of the internal Froude number 
(Fj = U/(g A p/p H) 54 ) less than 1/ir, a disturbance, was propagated upstream from a barrier 
moving through a fluid of depth H. He noted that for this case his analytical flow predic- 
tions became invalid, since they were dependent upon the assumption of undisturbed 
upstream flow. 

More recent experiments have been performed by Pao (1967, 1968), Martin and Long 
(1968), Debler and Daily (1971), Browand and Winant (1972), and Laws and Stevenson 
(1972). The common features of these experiments are the use of a linearly stratified salt 
solution for the working fluid, and use of very slow obstacle speeds. All these investigators 
have reported upstream influence with a wavy velocity profile, thus confirming the early 
observations of Long (1955), who reported an upstream velocity profile consisting of 
alternate jets and stagnation regions. 

Theoretical 

The problem of stratified flow about an obstacle has also been studied analytically by 
several investigators. These studies generally have taken one of three theoretical approaches. 

The first approach assumes inviscid steady flow everywhere and a uniform velocity profile 
upstream, includes nonlinear effects, and is generally based on the use of Long’s equation; 
for example, Long (1953). For internal Froude numbers greater than l/7r, good agreement 
is obtained between theoretical and experimental results. For lower internal Froude num- 
bers, the experiment shows upstream blocking, a condition which violates the theoretical 
assumption of undisturbed upstream flow. It is therefore not surprising that this approach 
fails to predict the upstream flow conditions for F. < l/7r. 

By using the formulation of an inviscid, linearized initial value problem, the second approach 
permits study of the formation of the blocking column. However, introduction of viscosity 
into this formulation requires the assumption of a Schmidt number (Sc = v/D) equal to 
unity, which is far from accurate for the usual experimental conditions involving salt solu- 
tions for which Sc * 800. This approach is the one used by Bretherton (1967). 

The third approach is a steady viscous formulation for very slow, low Reynolds number 
flows. The fluid is assumed to be nondiffusive, and the low velocity assumption allows 
the nonlinear inertial terms to be neglected in the flow equations. This formulation is thus 
a balance of the gravitational and viscous forces. This approach was used by Pao (1967, 
1968), who used two-parameter perturbation expansions to solve for the flow field in front 
of and above a horizontal flat plate. He obtained good agreement with his experiment. 
Martin and Long (1968) and Brown (1968) have also arrived at solutions that predict 
upstream disturbance for slow stratified flow past a horizontal flat plate. Graebel (1969) 
solved for the flow field about a vertical flat plate by means of matched-asymptotic-expan- 
sion theory. His solution also predicts the existence of an upstream influence. Janowitz 
(1971) used an integral transform approach to examine the low Reynold’s number flow of 
a linearly stratified fluid about a vertical flat plate. His solution predicts the existence of 
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an upstream blocking column within which the fluid recirculates in closed streamlines. The 
experimental work of Browand and Winant (1972) appears to confirm this result. 

In a separate class of theoretical papers, Benjamin (1970), Keady (1971), and McIntyre 
(1972) have examined the phenomenon of upstream influence from a more abstract view- 
point to determine whether Long’s hypothesis of no disturbance far upstream is ever 
strictly correct. Benjamin concludes that an upstream disturbance exists in the form of a 
uniform, long wave extending to steadily increasing distances ahead of the body whenever 
lee waves are present downstream. Keady extended Benjamin Vwork to the case of inter- 
facial waves in a two-fluid system. McIntyre, whose paper is still the subject of considerable 
controversy, discusses a kind of upstream influence different from that considered here. 

His theory predicts an upstream influence related to the formation of lee waves and mani- 
fested by an increase in upstream velocity with respect to the body rather than a decrease. 

Geophysical Implications 

The flow of a stratified ocean or atmosphere over a barrier such as a mountain ridge usually 
will produce lee waves downstream and blocking upstream. The prevailing west wind flow- 
ing over the San Gabriel Mountains east of Los Angeles, for example, produces blocking 
which makes it more difficult for the region’s pollutants to be dispersed. 

Background of This Study 

Previous studies have dealt with the slow flow about a two-dimensional body of either a 
linearly stratified fluid or of a two-layer fluid. The Unear-density model is an idealization 
which has proven to be very useful in studies of atmospheric and oceanic flows. In the real 
atmosphere, however, the near-linear density gradient may be interrupted by one or more 
inversions; and in the real ocean, a thermocline of one or more sheets is usually present. At 
these locations, the density undergoes a sudden but nevertheless continuous change, for 
which neither the Unear nor the two-layer model is a valid representation. The real condi- 
tion is somewhere between these two extremes. 

Long referred to this situation in Part III of his well-known trilogy on stratified flows 
(1955). He suggested that for strong inversions the shock-wave type of phenomena 
(hydraulic jumps and drops), such as exhibited in the flow of a two-layer fluid over a 
barrier (Long 1954), would dominate. For weak inversions, he proposed that the flow 
would behave Uke the Unear model. He concluded, however, that the intermediate case 
would be highly nonhnear and therefore difficult to solve analytically. 

Davis (1967), while doing research for his doctoral dissertation on the stability of oscilla- 
tory internal waves, used an experimental model consisting of a layer of fresh water on top 
of a layer of salt water, with a thin transition region at the interface. He recognized, within 
the transition layer, the propagation of a new type of solitary wave which had not been 
reported in the literature. This type of solitary wave was new, in that it could be propa- 
gated in a fluid of infinite depth, whereas internal solitary waves were known previously 
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only as shallow-water phenomena. Davis added an investigation of this new wave to his 
dissertation (see also Davis and Acrivos 1967). Independently and at about the same time, 
Benjamin completed a theoretical study in which he predicted the existence of such waves 
(1967). 

Davis’s experience suggests that the flow behavior of a stratified fluid in which the density 
changes suddenly over a short distance cannot necessarily be inferred from previous work 
on linearly stratified or two-layer fluids. In view of the importance of sudden density 
changes in geophysical fluids (atmospheric inversions and ocean thermoclines), it was 
decided to undertake a study of flow about an obstacle by a fluid with a sudden but 
continuous density change. 

Outline of This Study 

Experimental 

A laboratory flow model consisting of a layer of fresh water on top of a layer of salt water 
with a thin transition region at~the interface was used for this study. A vertical flat plate 
of 2.54-cm height was towed through the interfacial layer at various speeds and for various 
density profiles. Upstream and downstream velocity profiles were measured by a flow 
visualization technique. A conductivity probe was used to determine the density profile 
before and after each run. In addition to the shape of the density profile, the present work 
differs from that of most previous experiments in another respect. Those workers whose 
theoretical work depended upon neglecting the inertial terms in the momentum equation 
considered only very slow flow speeds in their experiments. The top speed considered by 
Laws and Stevenson (1972), for example, was only 0.07 mm/s, corresponding to a Reynolds 
number of 0(1) and a Richardson number (Ri = -g/p dp/dz/(dU/dz) 2 ) of 0(1 0 4 ). For the 
present study, flow speeds as high as 4 cm/s were examined. The corresponding dimension- 
less parameters were Re = 0( 1 0 3 ) and Ri = 0( 1 ). 

Theoretical 

The flow conditions described above were formulated as a nonlinear initial-value problem 
in a viscous diffusive medium. The Boussinesq approximation, which neglects the effect of 
density variation on the nonlinear inertial terms, was not used. The initial density profile 
was expressed mathematically by fitting experimental measurements to the equation 
p = p (1 - cu tanh «z), where dJ is a dimensionless constant and of’ is a measure of the 
thickness of the density-gradient region. Formulation as an initial-value problem permitted 
theoretical study of the development of the upstream influence. The formulation is 
analytically intractable, and the problem was therefore solved numerically. Comparison of 
theoretical upstream velocity profiles with those measured experimentally shows good 
agreement, but the theory fails to predict the existence of the reverse lee jet which was 
observed downstream in the density-gradient layer. 


5 



THEORETICAL STUDY 


Formulation 


Dimensionless-flow equations will now be developed for the case of a vertical flat plate of 
height 2b moving to the right with velocity U through the density-gradient layer of a vis- 
cous, diffusive, stably stratified fluid whose initial density profile is of the form p - p 
( 1 - co tanh ocz). The flow is viewed from a Cartesian coordinate system attached to the 
plate (see figure 1), and the plate is of infinite extent in the y direction. 



Figure 1. Geometry of the Flow Field 


The continuity equation for a diffusive medium is 


bp 9 3 , 

— + — (up) + — (wp) = p • 
3t 3x 3z 


The incompressibility equation is 


( 1 ) 


3u 3w 
— + — = 0 . 
3x 3z 


( 2 ) 


Since this equation is not quite correct for the case of a diffusive fluid such as salt water 
(see Phillips 1969, p. 15), it introduces a slight error to the formulation. The momentum 
equations are 


where 


Du 

0p 

P Dt 

3x 

Dw 

0p 

P Dt 

3z 


+ pV 2 u 
- pg + pV 2 w 


D 3 3 3 

— = — + u — + w — 
Dt 3t 3x 3z 



3x 2 


+ 


j ) 2 

3z 2 


(3a) 

(3b) 
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The stream function \p is defined by 

dip dip 

u = — , w = - — 

3z 3x 

and the y component of the vorticity <3 = V X q” = (£, 17, f) is 

3u 3w 
3z 3x 

so that by combining (4) and (5a) 

V 2 P = V- 


(4) 


(5a) 


(5b) 


Subtracting the x derivative of (3b) from the z derivative of (3a) yields 


3p Du 
3z Dt 


3 p Dw 
3x Dt 


+ P 


3r? 3 9 

— + — (ui?) + — (wn) 
3t 3x 3z 


dp 2 
= g— ■ + pV v . 
3x 


( 6 ) 


The Boussinesq approximation, which will not be made here, consists of neglecting the first 
two terms on the left side of (6), and replacing the factor p in the third term by the average 
density p . 

The boundary conditions for the system (1), (4), (5b), and (6) for t > 0 are 

as lx I, Izl -*■ p -»■ p (1 - w tanh az), ip ->• - Uz 

dp 

and at x = 0, IzKb, — = 0 
3x 


(7a) 

(7b) 


dip 

dx 


dip 

= 0 , — = 0 . 
3z 


(7c) 


Since u, w, and rj are defined in terms of \ p by (4) and (5b), no independent boundary con- 
ditions are possible for these quantities. Although (4), (5b), and (6) can be easily combined 
into one unwieldy equation, this is not done here, to retain clarity and manipulative ease. 
The initial conditions are 


and 


p = p (1 - w tanh az) everywhere at t = 0 


(8a) 


dip dip 
dx dz 


everywhere at t = 0 . 


(8b) 


With the use of the density scale-length a" 1 as a characteristic length, U as a characteristic 
velocity, and p as a characteristic density, these equations can be put into dimensionless 
form. Thus 
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( 9 ) 


<x,z) = « (x, z), (u, w) = — (u, w), p = 3 

U p 


Uat, 


V = — V, 
Ua 


~ Q 

* “TT* 


Equation (1) then becomes 


where 


Equation (6) becomes 


3p a 3 1 

— + — (up).+ — (wp) = V P 


3 1 3x 


3z 


Sc • Re 


Sc, Schmidt number = vjD 
Re, Reynolds number = (U/cu>) . 


(10) 


3/d DiT dp Dw ^ 

— — ' — — + P 

3z DT 3x DT 


3p 3 ^ ^ 3 ^ ^ 

— + — (up) + — (wr?) 
37 3x 37 


1 3p 
F 2 3x 


+ —V P, 
Re 


where 


F, Froude number = V/fg/a)^ . 


01 ) 


Similarly, (4) and (5b) become 


and 


~ 3^ -9^ 

u = — , w = 

37 9x 


V 2 \p = V ■ 


( 12 ) 


(13) 


Note that the Froude number appearing in (1 1) is of the overall form, and is not the inter- 
nal Froude number, F. = U/(g A p/p H) % , that commonly appears in the formulation of 
stratified flow problems. This is a result of treating density in the present time-dependent 
problem as a separate unknown with a separate governing equation rather than com- 
bining the density profile into the momentum equation. To facilitate comparison of pre- 
vious work with the results of this study, internal Froude numbers were computed for the 
conditions of the experimental runs. 

The dimensionless form of the boundary conditions for t >0 is 


as lx" I , I7I-*- 00 , p' ->■ (1 - w tanh7) , ip -+- 7 


(14a) 


9 



and at x = 0 , 


(14b) 


Iz I < ab , 




dx 


dip 

— = 0. 
dz 


(14c) 


The initial conditions become 

p = (1 -co tanh z) everywhere at T = 0 

and 

dip dip „ , 

— = — - 0 everywhere at t = 0 . 

dx dz 

Henceforth, the tildes (~) will be omitted and all variables will be understood to be 
dimensionless. 

Numerical Solution 

Algorithm for Solution 

The solution of a modified form of (10) through (15) was accomplished by writing them as 
finite difference equations and solving these on a digital computer according to the following 
algorithm: 

(a) At t = 0, the initial conditions (15a) and (15b) exist. 

(b) At t = 0 + , the flow is irrotational, with t? = 0. The solution for the \p field is 
accomplished by making an initial guess of ip = -z everywhere and then conver- 
ging to the correct result by overrelaxation of (13), subject to the boundary 
conditions (14a) and (14c). 

(c) The u and w velocity fields at t = 0 + are computed from (12). 

(d) Equation (10), subject to the boundary conditions (14a) and (14b), is used to 
make a time step in the density field. 

(e) Equation (1 1) is used to make a time step in the vorticity field. [The boundary 
conditions on ip result in the following conditions on tj: as lx I, Izl ->• °°, 17 -> 0; 
and at x = 0, Izl < ab, (u = d\p/dz = 0) => 7? = 3 2 i///9x 2 ] . 

(0 With the new result for the vorticity field, a new value for the ip field is computed 
by overrelaxing (13), subject to (14a) and (14c). 

(g) The u and w velocity fields are computed from (12). Steps d, e, f, and g are 
repeated until a specified time t = t f has been reached. The computer can be 


(15a) 


(15b) 
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requested to print out the p, \p, rj, u, and w fields at t = 0 + , t = t f , and at inter- 
mediate times, so that the time development of the flow field can be studied. 

As mentioned above, the numerical solution is accomplished with a modified form of 
(10) through (15). The nature of these modifications will now be discussed. 

Modification of Nonlinear Boussinesq Terms 


Numerical evaluation of nonlinear terms was accomplished by a technique called the special 
three-point noncentral differences method (Torrance and Rockett 1969). To use this 
method, it is necessary that nonlinear terms be of the form 3 (fg)/dx rather than of the 
form g (9f/9x). This is accomplished in the case of the first two terms on the left side of 
(1 1) by recalling (2) and adding u (3u/3x + 3w/3z) to the factor Du/Dt and w (9u/3x 
+ 3w/3z) to the factor Dw/Dt. Equation (11) then becomes 


dp 

dz 


+ P 


dp 

3w a a 

) - — 

— + — (uw) + — (w 2 ) 

J 3x 

3t dx dz 


f3u 9 d 1 dp 

rar, a a 1 i dp i 

— + — (UT?) + — (wi?) = — — + —V v - 

[_at 9x 3z J p 2 3x Re 


Stretching the Coordinate System 


(16) 


Previous attempts at numerical solution of similar problems have had difficulty with the 
application of boundary conditions at infinity. A dilemma is posed by the conflicting 
requirements of keeping the finite-difference grid mesh fine in the vicinity of the body 
where the gradients are largest, and at the same time having the last grid point be far 
enough away from the body that it is an accurate approximation of infinity. One is forced 
to choose among an overly coarse mesh near the body, or prohibitively large computer 
core requirements, or an infinity so close to the body that errors are-introduced. This 
problem was circumvented by writing (10) and (1 2) through (16) with respect to a 
coordinate system such that its mapping to real-world coordinates has a fine-grid mesh 
near the origin, a coarser mesh farther away, and a last grid point as close to infinity as 
desired. This mapping, which involves the hyperbolic tangent function, will now be 
developed. 


The coordinate system sought (x, z) is one whose mapping to ordinary Cartesian coordinates 
(x, z) is such that as x" and z vary from -1 to +1 , x and z vary from to +°°. The trans- 
formation for such a system is given by 

x = tanh ax, z = tanh cz (17) 


where 


— = a • sech 2 ax = a(l - tanh 2 ax) = a(l - x 2 ) (18a) 

dx 
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and 


dz 

-= c(l 
dz 


-2 




Therefore, for any arbitrary function f, 


and 


3f dx df 3f 

— = = a (1 - x 2 )- 

3x dx 3x 3x 


3f , 3f 

= C (1 - Z 2 ) — 

3z 3z ’ 


so that from (12) 


(18b) 


(19a) 


(19b) 


and 


dd/ „ 3 \b 

u = — = c (1 - z 2 ) — = c (1 - z 2 ) u 
3z 3z 


(19c) 


Finally 


3i// di// _ _ 

w = = -a (1 -x 2 ) — = a (1 -x 2 ) w. 

3x 3x 


(19d) 


and 


3 2 f _ 3 

3x 2 3x 


fd - x 2 ) -1 = a 2 

L dx. 


• 2 

3x 2 


a‘(l -x 2 ) — - 2a 2 x (1 -x 2 )— ( 1 9e) 


3f 

3x 


3z 2 dz 2 dz 


(190 


At this point, it is useful to restate the untransformed dimensionless governing equations. 


3p 3 3 

— + — (up) + — (wp) = 
3t 3x 3z 


dp 

dz 


3^ 

3z ’ V 

V 2 <p = V- 


- — T v2 p 
Sc • Re 

(20a) 

> Tdw 3 3 ."I 

- — + — W + — (w : 2 ) 

: [3t 3x 3z ] 


1 dp 1 

— — + — V 2 v 
p 2 3x Re 

(20b) 

1 

QP 1 CP 
X 1^. 

(20c) 

V ■ 

(20d) 
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Boundary conditions for t > 0 are 


as lx I, Izl -»■ °o, p -> (1 - cu tanh z), \p ->■ -z 


and at x = 0, Izl < ab, — = 0 

9x 


d\ p 9i p 

— = 0 , — = 0 . 

9x 9z 


Initial conditions are 


p = (1 - oj tanh z) everywhere at t = 0 


d\p 9i p 

— = — = 0 everywhere at t = 0 . 

9x 9z 


(20e) 

(2 Of) 

(20g) 

(20h) 

(20i) 


Substitution of the transformation (19) into the governing equations (20) results in a sys- 
tem of equations which can be put into finite-difference form and solved numerically by 
the algorithm discussed above. The transformed equations are so lengthy that their presen- 
tation has been deferred to appendix 1. The finite-difference form of these equations is 
given in appendix 2, and the computer program is presented in appendix 3. 

At this point, it is appropriate to mention another approximation necessitated by numeri- 
cal modeling. Although the vertical flat plate used in the experimental study was very thin 
in the x-direction compared to its height, it was necessary that the numerical model of the 
plate be at least two x-grid widths thick to avoid instability in the solution caused by singu- 
larities at (0, b) and (0, -b). This approximation is not believed to have significantly 
affected theoretical results upstream, but was found to be responsible for the discrepancy 
between theoretical and experimental results downstream. 


Numerical Results 

Figure 2 presents a typical theoretical result for the velocity profiles at four upstream 
locations at some time, t = tj . Figure 3 shows a typical theoretical prediction for the 
variation of center line velocity with upstream distance at t = t x . Further information on 
numerical results will be presented in the section, Discussion of Results, where they are 
interpreted in the light of experimental findings. 
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EXPERIMENTAL STUDY 


Description of Experiments 

Filling the Channel 

The experiments were carried out in a Plexiglass channel 10 m long, 60 cm deep, and 
35.5 cm wide. On the day before a set of experimental runs, the channel would be filled 
to a depth of 23 cm with a uniform salt water solution lightly tinted with nigrosine. To 
ensure its uniformity, the solution in the channel was thoroughly agitated and then allowed 
to sit for at least 12 hours. Two salt solution densities were used during this study, 

1.013 g/cm 3 and 1.026 g/cm 3 . On the day of the experimental runs, a layer of fresh water 
was slowly floated onto the salt solution until the total fluid depth was 46 cm. This was 
accomplished by letting the fresh water from the filling hose flow into a sponge rubber 
basin which was floating on a Styrofoam raft. The water seeped through the sponge, 
thereby dissipated most of its momentum, and floated with little mixing onto the salt 
solution. Gradient-layer thicknesses of the order of 4 cm were achieved by this method. 

By making certain that the temperature of the inflowing fresh water differed by no more 
than 2° C from that of the salt water onto which it was floated, convective fluid motions 
were reduced. These small secondary motions could not be entirely eliminated, however. 

Measurement of Density Profiles 

Density profiles were measured before and after each run with a conductivity probe. The 
probe was calibrated with standard salt solutions whose densities had been measured by a 
hydrometer to an accuracy of ±0.0005 g/cm 3 . Density measurements made with the probe 
are estimated to be accurate to ±0.0005 g/cm 3 , while the z elevations for these density 
measurements have an estimated accuracy of ±0.2 cm. The resulting density profiles were 
fitted to Benjamin’s expression ( 1 967), p = p ( 1 - 63' tanh az), where p = 1/2 (p u + P e ) 

= average density of upper and lower layers; dJ = (p e - p u )/(p e + p u ); a is a number deter- 
mined by least-squares fit whose reciprocal is a measure of the gradient thickness, and z is 
the height above the level at which p = p . Figure 4 shows a typical experimentally-meas- 
ured density profile and the function to which it was fitted. Interestingly enough, it was 
found that the increase in the gradient layer thickness a' 1 was very nearly a linear function 
of the number of transits made through the layer by the vertical flat plate. This relation- 
ship is illustrated in figure 5. Since the diffusivity of salt in water is very low, the thickening 
is almost entirely due to mixing. 

Typically, about four experiment runs and four subsequent retrievals of the plate could be 
made during the course of a day. Approximately an hour was allowed between runs for 
residual fluid motions to dampen out. 
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Figure 5. Thickening of the Gradient Layer Caused by 
Transits of the Vertical Flat Plate 
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Plate and Towing System 

The vertical flat plate was a Plexiglass strip 2.54 cm high, 35.0 cm wide, and 2 mm thick. 
This was suspended below a carriage by two Plexiglass struts which had been machined to 
an easy shape to minimize their effect on the flow. The carriage rolled on four double- 
race ball-bearing roller skate wheels along rails made from aluminum angle. The towing 
mechanism consisted of a variable-speed motor at one end of the channel and a system of 
counterweights at the other end. The motor and the counterweights were connected to the 
carriage over pulleys by size 20 Hygrade flexible fiberglass cord.* Of several types of cords 
and cables tried, this material proved to be the most satisfactory for minimizing jerkiness 
in the towing motion. Considerable time was spent in trial and error with the towing sys- 
tem in an attempt to achieve smooth towing motion. Although towing motion was made 
adequately smooth for the range of plate speeds considered in this study, it was less than 
adequate at very slow speeds. During an experimental run, the motor was reversed so that 
the cord unwound from the motor pulley and the counterweight pulled the carriage along 
the rails. After each run, the plate was retrieved to its original starting position in prepara- 
tion for the next run. Figure 6 is a schematic drawing of the plate and towing system. 



Figure 6. Plate and Towing System 


Early in the study, some qualitative experiments were run with two other body shapes: 
a 2.54 cm-diameter circular cylinder, and a 2.54 cm-thick lenticular cylinder. The vertical 
flat plate was chosen for all quantitative experiments, because it produced the strongest 
upstream influence. 


* The cord is manufactured by Maikel & Sons of Norristown, Pa. 
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Flow Visualization 


Visualization of fluid motion in the gradient layer was aided by injecting droplets of a 
mixture of carbon tetrachloride, mineral oil, and oil red dye into the upper layer with a 
syringe and hypodermic needle. The density of the mixture was such that the droplets 
were neutrally buoyant at the middle of the gradient layer. The syringe and hypodermic 
needle were required both to control droplet size and to prevent the mixture from being 
held at the free surface by surface tension. Droplets of a mixture containing toluene 
instead of mineral oil were tried earlier in the study. That mixture was rejected because 
of the toxicity of toluene fumes, high volatility of the toluene which made it difficult to 
control the mixture density, and the mixture’s tendency to dissolve and permanently 
discolor the Plexiglass channel. 

In addition to the droplets in the gradient layer, upstream flow visualization was accom- 
plished by the use of vertical dye lines. These dye lines were created by dropping 
nigrosine dye crystals upstream of the approaching body. Arrival of the upstream distur- 
bance was made visible by the horizontal displacement of the dye lines. When the first 
set of dye lines became so distorted that accurate measurement of their displacement was 
no longer possible, a second group of dye crystals was dropped. The second drop created 
dye lines which made it possible to evaluate the flow in the near upstream vicinity of the 
plate. Figure 7 is a photograph showing the displacement of dye lines due to the arrival of 
an upstream disturbance. 

Measurement of Upstream Velocity Profiles 

Upstream velocity profiles were evaluated by photographing the droplets and dye lines at 
known time intervals and measuring the displacement which they experienced during those 
intervals. As figure 7 shows, the displacements were measured with respect to a rectangular 
grid fixed to the rear wall of the channel. Also visible in the figure is one of the three 
stopclocks used. These clocks were started simultaneously with the start of the plate 
towing motor so that they measured the total duration of the run as well as the time 
intervals between photographs. The clock shown, an Aristo stopclock, measures time in 
minutes and seconds. Also used, but not shown in the figure, was a larger stopclock 
(manufactured by the A. W. Haydon Co.) which measures time in seconds and hundredths 
of seconds. The time intervals between photographs were measured with an accuracy of 
±0.01 seconds. 

As shown schematically in figure 8, two cameras with overlapping fields of view were used 
so that detailed photographs were made of a section of the channel about 125 cm long. 
Nikon-F cameras were used with 135 mm Nikkor telephoto lenses. Each camera had a 
motorized shutter release and film advance which made possible short time intervals 
between successive exposures. Camera motors were connected in parallel to a single switch 
so that the cameras made simultaneous exposures. The cameras were focused on the mid- 
plane of the channel, so that both the rectangular grid and the clocks were in focus. 
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Figure 8. Placement of the Cameras 

Typical shutter and aperture settings were 1/250 second at//4.0, using Hi-Speed Ekta- 
chrome color film developed for an effective film speed of ASA 400. 

Measurement of Downstream Velocity Profiles 

The measurement of downstream velocity profiles was considerably less accurate, not only 
because no dye crystals were dropped there, but also because the mixing in the wake com- 
bined with the sudden change of index of refraction in the gradient region to cause such 
optical distortion that accurate displacement measurements were impossible. By observing 
the motion of the oil droplets, however, qualitative assessments of the downstream flow 
could be made. The most striking feature of the downstream flow was the presence of a 
reverse jet in the gradient layer, which carried oil droplets about 75 cm downstream with 
respect to the fixed cameras before dissipating. This reverse jet is discussed in the section, 
Discussion of Results. 

Measurement of Body Speed 

The speed of the vertical fiat plate was determined by two methods. By stopping the 
towing motor and the clocks simultaneously at the completion of a run, the body speed 
could be computed directly from the length and duration of the run. Secondly, the body 
speed could be computed as the slope of a straight line determined from a least-squares fit 
of photographically recorded body positions and elapsed times. Excellent agreement was 
obtained between the two methods. 
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Measurement of Speed of Disturbance Front 

The least-squares method was also used to determine the propagation speed of the front of 
the upstream disturbance. 

Parameters for Experimental Runs 

Table 1 lists values of the density and flow parameters for the ten quantitative experimen- 
tal runs. Data from the last six runs are considered most reliable. 



Flow Parameters for Experimental Runs 
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DISCUSSION OF RESULTS 


Qualitative Description of Observed Flow 


Figure 9 is a schematic drawing which shows the type of flow observed to be caused by the 
motion of a vertical flat plate in a region of sudden density variation. 







' SOLITARY 
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REVERSE LEE JET 
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UPSTREAM INFLUENCE 
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Figure 9. Schematic Drawing of the Type of Flow Phenomena Induced by the Motion of a 
Vertical Flat Plate in a Region of Sudden Density Variation 

At the front of the upstream disturbance an internal solitary wave of the type predicted by 
Benjamin (1967) and observed by Davis (1967) was always present. Behind the solitary 
wave was a region of upstream influence. This upstream influence did not take the form of 
complete blocking such as has been observed in the very slow motion of bodies through 
stratified fluids. Nevertheless, a strong upstream wake was observed in the gradient layer. 
(The solitary wave and the upstream influence when discussed together are called here the 
upstream disturbance.) Both the solitary wave and the upstream influence propagated up- 
stream to steadily increasing distances in front of the plate. This is consistent with the 
internal wave concept of upstream influence, discussed in the Introduction. The connection, 
if any, between the upstream influence and the solitary wave is not entirely clear. Present 
indications are, however, that the solitary wave is a distinct phenomenon generated by the 
impulsive start of plate motion. 

As shown also in figure 9,, the upstream flow in the fresh- and salt-water layers was opposite 
in direction to that in the gradient layer. 

The downstream flow observed in the gradient layer was somewhat surprising. Nearby the 
body, the flow looked like ordinary wake flow. About 10 cm downstream, however, a 
reverse jet was observed to carry oil droplets downstream not only with respect to the 
moving plate, but also with respect to the fixed cameras. The following explanation is 
suggested: The observed upstream influence did not take the form of complete blocking, 
and some upstream fluid was sufficiently displaced vertically that it cleared the plate. 

After passing the plate, the stratified fluid rushed back to its position of neutral buoyancy. 
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In doing so, it accumulated a residual rearward momentum manifested in the form of the 
observed lee jet. 

The reverse lee jet dissipated about a meter downstream where it was replaced by a system 
of lee waves in which the oil particles, as seen by the fixed cameras, simply orbited in 
closed paths. 


Examination of the Internal Solitary Wave 

Figure 10 is a photograph showing the propagation of an internal solitary wave from left to 
right ahead of the upstream influence which is in turn ahead of the vertical flat plate. 


Comparison With Theory 

The solitary waves which were observed in this study to head the upstream disturbance 
have been examined in the light of the work of Benjamin (1967) and Davis and Acrivos 
(1967). 

Benjamin predicted that for a region whose density profile is described by p = p (1 - c3 
tanh az), an infinitesimal solitary wave of the first mode would propagate with the speed 
c° w = (goJ/2a) Vi . A wave with finite amplitude a is predicted to propagate with the speed 


3 

(1 +-ota) 


'h 


( 21 ) 


In the present study, the speed c^ could be determined accurately from photographs by 
the least-squares method (see the section, Experimental Study). Thus, it was possible to 
compare Benjamin’s theory with experimental results by measuring the observed solitary 
wave amplitudes, using (21) to predict wave speeds, and comparing these predictions with 
experiment. For the purpose of these measurements, the amplitude was defined as the 
maximum vertical displacement of the streamline at z = a' 1 . The results of the comparison 
are presented in table 2. 

Considering the crudeness of the experimental amplitude measurements, the agreement ' 
between the last two columns of table 2 is quite good. 

A correlation between the dimensionless amplitude, aa, and the dimensionless group 
X = g/2ac^, £n (p s /p u ) is predicted by Davis’s theory. Moreover, for small values of aa, 
the relation X = 2.0 - 1 .2 aa is predicted to be a good approximation to this correlation. 
Figure 1 1 is a plot of X vs aa for the solitary waves observed in the present study. The 
theoretical predictions of both Davis and Benjamin are plotted in the same figure. The 
agreement between these theories and the results of the present experiment is quite good. 
Since the theories have nothing to do with upstream influence, this agreement indicates 
that the observed solitary waves are generated by the impulsive start of plate motion and 
are independent of the existence of upstream influence. 
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Figure 10. Propagation of an Internal Solitary Wave from Left to Right Ahead of the Vertical Flat Plate 



: Table 2 


Comparison of Benjamin’s Predictions for 
Internal Solitary-Wave Speed With 
Experimental Measurements 


Run 

CO 

a 

(cm' 1 ) 

c o _/W Y /2 

sw \2 a) 
(cm/s) 

Measured 

Amplitude 

a 

(cm) 

Measured 
Wave Speed 
c 

S W 

(cm/s) 

Calculated 

Wave Speed 
c 

SW 

[Equation (21)] 
(cm/sec) 

13 

0.00771 

0.50 

2.75 

2.5 

3.73 

3.64 

14 

0.00771 

0.48 

2.81 

0.6 

2.95 

3.04 



2.44 

1.69 

1.5 

3.18 

3.02 

17 

0.01417 

1.09 

2.52 

2.8 

4.75 

4.24 

19 

0.01417 

0.73 

3.08 

1.8 

4.11 

4.12 

20 

0.01417 

0.65 

3.27 

1.3 

4.06 

4.01 

21 

0.01408 

0.32 

4.64 

0.0 

4.62 

4.64 

22 

0.01408 

0.31 

4.72 

0.0 

4.65 

4.72 

23 

0.01408 

0.30 

4.80 

0.6 

5.02 

5.05 

24 

0.01408 

0.29 

4.88 

0.6 

5.05 

5.13 
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Correlation of Plate Velocity and Wave Speed 

A plot of solitary wave speed c versus plate speed U appears in figure 12. The correlation 
between these two quantities is close enough to suggest a contradiction with the above 
conclusion that the solitary wave is independent of the existence of upstream influence. 
This writer suggests, however, that the correlation shown in figure 12 is quite indirect: 
that wave speed c m depends on the wave amplitude a, which depends on the strength of 
the starting impulse of the plate and on the speed U for which the towing motor is set when 
it is switched on. 

The results are therefore not as clear as one might wish, but the bulk of the evidence seems 
to indicate that the initially generated solitary wave and the upstream influence do not 
interact. 

Examination of Upstream Influence 
Brunt-Vaisala Frequency 

Now, recall from the Introduction the significance to the propagation of internal waves of 
the Brunt-Vaisala frequency N = (-g/p dp/dz) Vl . It was noted that: 

(a) For stably stratified fluids, dp/dz < 0, so that N>0. 

(b) No internal wave can be propagated with frequency higher than the local value 
of N. 

(c) The direction of propagation of a wave of frequency a is cos' 1 (a/N) from the 
vertical. 

Finally, it was pointed out that one may consider upstream influence to be the propagation 
of an internal wave having zero frequency, and that the direction of propagation of such a 
wave would be horizontal. 

Figure 13 is a plot of N versus elevation az for a typical density profile of the type consi- 
dered in this study. Note that N reaches its maximum at the center line az = 0 and drops 
off rapidly to zero at short distances above and below the centerline. The significance of 
this is that any wave having a frequency o such that 0 < a < N would propagate in an 
increasingly vertical direction as it left the centerline, but could not propagate beyond the 
elevation at which o = N. Thus, such waves would be trapped within the gradient layer 
and after multiple reflections would eventually die out. The waves with higher frequencies 
would be reflected more often and would tend to die out faster. 

Phase Velocity of the Upstream Influence 

Experimental measurement of the propagation speed c^ of upstream influence was not 
possible in this study due to the presence of the solitary wave. It was not clear where the 
disturbance caused by the solitary wave ended and where the front of the upstream 
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Figure 13. Brunt-Vaisala Frequency N 
as a Function of Elevation above 

Centerline in a Region of Sudden Figure 14. Correlation of Predicted Upstream- influe nee Propagation Speed c u . with 

Density Variation Froude Number F 








influence began. Moreover, the velocity profiles were not observed over a time interval long 
enough to permit experimental evaluation of their variation with time. A phase velocity for 
the upstream influence could be determined, however, from the results of the numerical 
study. Since, as will be shown later, the numerical solutions successfully predicted upstream 
velocity profiles as a function of time and position, the determination of c^ from these 
solutions is thought to be quite reliable. Table 3 presents these results. In figure 14, c^/U 
is plotted as a function of the Froude number F = U/(g/a) H . 

Suggested Relationship Between c u . and Group Velocity of Infinitesimal Internal Waves 

The i th component of the group velocity of an infinitesimal internal wave of frequency o 
and wave number vector k is (c g ). = 3o/3k., i = 1 , 2, 3. The relationship o - o (k.) 
between frequency and wave number is called the dispersion relation, and is obtained from 
the analytic solution of an eigenvalue problem formulated from the equations of motion 
for a stratified fluid with appropriate boundary conditions (see Phillips 1969, pp. 161-165). 
Phillips points out that for the general case in which the Brunt-Vaisala frequency is an 
arbitrary function of z, simple analytic solutions to the eigenvalue problem are not possi- 
ble. For the present study, N(z) = (gc3a) Vl sech a z (see figure 13), and the solution of the 
eigenvalue problem is such' a formidable task that it has not been attempted here. Never- 
theless, it is tempting to suggest that the speed c^ is related to the group velocity of 
infinitesimal internal waves. The existence of such a relationship is suggested by the results 
of Kao, Pao, and Wei (1972), who found that when a linearly stratified fluid in a channel 
suddenly begins to discharge into a line sink, a disturbance propagates upstream from the 
sink with velocity equal to the group velocity of internal gravity waves of zero horizontal 
wave number. In view of the analogy between rotating and density stratified flows, a 
relationship between c^ and group velocity is also suggested by the study of Benjamin and 
Barnard (1964). They found that the motion of a cavity in a rotating fluid induced an 
upstream disturbance which propagated at the maximum group velocity for infinitesimal 
traveling waves. Confirmation of the existence of such a relationship for the conditions of 
this study will await the solution for the dispersion relation. 

Comparison of Observed Upstream Influence With Predictions of Janowitz and Experi- 
ments of Browand and Winant 

The Introduction mentioned the theoretical study by Janowitz (1971) of the slow motion 
of a vertical flat plate of height 2b in a linearly stratified fluid. This study, whose results 
are valid when 1/Ri <K Re 1, predicts a steady-state, upstream blocking column of length 
10~ 2 b (Re • Ri) composed of cells of closed streamlines. The experiments of Browand and 
Winant (1972) appear to confirm the prediction of a region of recirculating flow upstream. 

In the present study, the body speeds were much higher, with Ri = 0(1) and Re » 1. Up- 
stream influence was definitely present, but not in the extreme form of complete blocking. 
Nevertheless, it is instructive to calculate the slug lengths 10‘ 2 b (Re • Ri) for the conditions 
observed. These are presented in table 4. 
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Table 3 


Propagation Speed of Upstream Influence 


Run 

Body 

Speed 

U 

(cm/s) 

Froude 

Number 

F 

U/(g/a)* 

Solitary 

Wave 

Speed 

c 

S W 

(cm/s) 

Propagation 
Speed of 
Upstream 
Influence 

c . 

U 1 

(cm/s) 

c ./U 

13 

3.18 

0.0718 

3.73 

0.225 

0.07b 

14 

1.95 

0.0432 

2.95 

0.857 

0.440 

16 

1.94 

0.0968 

3.18 

0.058 

0.030 

17 

3.80 

0.1268 

4.75 

0.080 

0.021 

19 

2.60 

0.0710 

4.11 

0.216 

. 0.083 

20 

2.60 

0.0670 

4.06 

0.325 

0.125 

21 

3.31 

0.0598 

4.62 

0.724 

0.219 

22 

3.16 

0.0558 

4.65 

0.800 

0.253 

23 

4.08 

0.0709 

5.02 

0.354 

0.087 

24 

4.15 

0.0708 

5.05 

0.324 



0.078 
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Table 4 


Length of Upstream Blocking Column as Predicted by the 
Janowitz Formula for Present Experiments 


Run 

13 

14 

16 

17 

19 

20 

21 

22 

23 

24 

10~ 2 b (Re • Ri) 
(cm) 

3.0 

5.0 

48. 

11. 

11. 

10. 

4.0 

4.0 

3.0 

3.0 


These results confirm that Janowitz’s theory is indeed inapplicable to this study, for no 
recirculating flow of the type described was observed within the predicted distances of the 
plate. Had such closed streamline flow been present, the oil droplets and dye lines would 
surely have revealed it. 

Influence of the End Wall 

Another way in which this study differs from the experiment of Browand and Winant is 
that it was unnecessary to correct the data for the finite length of the channel. The front 
of the upstream disturbance was still several meters from the end wall when the velocity 
profile photographs were taken. 

Theoretical and Experimental Results for Upstream Flow Field 

The upstream velocity profiles predicted by the numerical study (see section, Theoretical 
Study) will now be compared with the experimental results. Figures 15 and 16 present the 
results for two typical experimental runs together with the corresponding theoretical 
results. Profiles of dimensionless velocity with respect to the plate, u/U, are shown at four 
upstream locations. Agreement between theory and experiment is excellent in the gradient 
layer. Theoretical flow velocities in the homogeneous layers are in good qualitative agree- 
ment with the experimental, although the predicted velocities are somewhat greater than 
those actually observed. 

Figures 17 and 18 show theoretical and experimental results for the variation of centerline 
(z = 0) velocity with upstream position, ax. Here again, the agreement between theory and 
experiment is excellent. Shown on the same plots are the predictions for upstream-velocity 
variation for potential flow. These plots show that the upstream velocity variation is 
somewhat like potential flow in the immediate vicinity of the plate, but rapidly deviates 
from this idealization for ax > 0.5. 

Finally, these experimental results are compared with the theoretical predictions of 
Janowitz in figure 19, where dimensionless centerline velocity with respect to the plate, 
u/U, is plotted versus Ix/L 3 V '* . The quantity L, which Janowitz defines as L = b/(Ur , /)3g) 1/3 , 
can be expressed in terms of the parameters used in the present study as follows. 
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Figure 17. Theoretical and Experimental Results for the Variation of Centerline Velocity 

with Upstream Position, Run 21 . 

Janowitz considered the density profile 
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This study has considered 
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Figure 18. Theoretical and Experimental Results for the Variation of Centerline Velocity 

with Upstream Position, Run 23 

Thus, by analogy between (22b) and (23c), one can define for a plate of height 2b 

'< L = b/(UiVc3ag) 1/3 (24) 

Figure 19 again illustrates that Janowitz’s theory, which is formulated for high Richardson 
number and very low Reynolds number, does not apply to the present experiment. Agree- 
ment is particularly poor near the body, where Janowitz’s theory predicts the existence of 
a blocking column with recirculating flow. 

Approach to Steady State 

In figures 20 and 21 are plotted, for the two runs and at the four upstream stations pre- 
viously considered, the variation of the theoretical centerline velocity u/U with dimension- 
less time Uat. The times during which the experimental data were taken are indicated in 
the figures. These plots indicate that at the time of the experimental measurements the 
flow near the body had nearly reached steady state, but that farther upstream it was still 
quite unsteady. 
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Figure 19. Comparison of Experimental Results for Positional Variation of 
Centerline Velocity with Janowitz's Theory for Slow-flow Blocking 



Figure 20. Approach to Steady State. Dimensionless Centerline Velocity as a 
Function of Dimensionless Time at Four Upstream Locations, Run 21 
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Figure 21. Approach to Steady State. Dimensionless Centerline Velocity as a 
Function of Dimensionless Time at Four Upstream Locations, Run 23 

Relative Importance of Various Terms in the Theoretical Formulation and Conclusions 
about the Nature of the Observed Upstream Influence 

It will be recalled that the formulation presented in the section, Theoretical Study, was 
made as general as possible. It treated the flow of a diffusive, viscous, non-Boussinesq, 
stratified fluid about a vertical flat plate. A series of numerical solutions were run on the 
computer to evaluate the relative importance of various terms in the governing equations 
and to determine thereby the nature of the observed upstream influence: 

(a) When the diffusive term 1 /(Sc • Re) v 2 p was neglected from the formulation, 
the numerical results for the velocity field were absolutely unchanged. This is 
not surprising, since, for the conditions of these experiments, Sc = 800, while 
Re varied from 84 to 1538. Thus, for the conditions of the present study, the 
flow was nondiffusive. 

(b) The Boussinesq approximation resulted in a small but insignificant (1 percent) 
decrease in the predicted velocities. Evidently, the Boussinesq approximation 
is a valid assumption for the flow under study. 

(c) Neglecting the viscous terms 1/Re V 2 r? also resulted in a negligibly small 

(1 percent) decrease in the predicted velocities. Hence, one concludes that the 
plate speeds were fast enough that the upstream flow could be regarded as 
in viscid. 
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(d) Variation in the Froude number F = U/fg/a) 14 , however, had an overwhelming 
impact on the numerical solution. An increase in F resulted in a dramatic 
reduction in the length of the predicted region of upstream influence. 

From these results one can conclude that the flow under study was governed almost 
entirely by a balance between inertial and gravitational forces. This is a significant result, 
for it means that the observed upstream influence is of an entirely different regime from 
the slow-flow blocking studied by Janowitz and others. In the case of these very slow 
flows, a balance exists between the viscous and gravitational forces. 

The effect of the Froude number on the numerical solution is illustrated in figures 22, 23, 
and 24. Figure 22 shows the predicted variation of centerline velocity with upstream posi- 
tion for F = 0.02, 0.06, 0. 10, and 0.20 and for potential flow. This figure shows that as F 
increases, the upstream influence becomes weaker and the variation of velocity with posi- 
tion approaches that for potential flow. For very small values of F, on the other hand, the 
upstream solution shows a closed-streamline recirculation region near the body, a result . 
suggestive of the findings of Janowitz. This causes the solution to be somewhat unstable, 
as is shown clearly in figure 24, where dimensionless velocity at four upstream positions is 
plotted versus dimensionless time for F = 0.02. This result contrasts with the very stable 
solution for F = 0. 10, as illustrated in figure 23. For intermediate values of F corres- 
ponding to the conditions of our experiments, the numerical solutions were quite stable 
with no evidence of recirculating flow near the body. 

Examination of Downstream Flow Field 

The photograph presented as figure 25 shows a dye line created by a crystal that happened 
to drop downstream of the moving plate. The shape of the dye line in the gradient layer 
shows very clearly the existence of the reverse lee jet previously described. 

As mentioned earlier, there were difficulties in obtaining downstream velocity profiles. The 
optical distortion in the downstream-gradient layer is evident in the figure. Nevertheless, 
an attempt was made to determine downstream-velocity profiles from the motion of the 
oil drops that were not swept upstream by the plate. This attempt is summarized in figures 
26 and 2.7. Figure 26 presents profiles of dimensionless velocity with respect to the plate 
u/U at four downstream locations. Figure 27 shows the variation of centerline velocity 
with downstream position ax. 

The numerical solution, which models a plate of finite thickness, failed to predict the 
existence of the reverse lee jet. When the numerical model of the plate was shrunk to zero 
thickness in the x-direction, however, the solution for Uat < 40 did indeed show the begin- 
ning of the development of a reverse downstream wake. Unfortunately, the singularities at 
the tips of the infinitesimally thin plate created a disturbance which caused the numerical 
solution to become noticeably unstable both upstream and downstream for Uat = 40 and 
to become increasingly degenerate with time. The problem of the numerical solution of 
the equations of motion near mathematically sharp edges has been studied by Lugt and 
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Figure 22. Effect of Froude Number on the Numerical Solution, Showing Variation 
of Centerline Velocity with Upstream Position 


Ohring (1971). They found that they could avoid numerical instability only by manipu- 
lating their coordinate system so that the singularity did not coincide with a grid point. 
This is not a feasible alternative for the geometry of the present study. 

Thus, although qualitative prediction of the nature of the downstream flow could be 
accomplished by numerical modeling of a plate of infinitesimal thickness, that procedure 
created more problems than it solved. 

An Examination of the Numerical Study 

In summary, the numerical model successfully predicted the observed flow, including: 

(a) The wavelike propagation of the upstream influence to steadily increasing 
distances in front of the plate 

(b) The correct shape and magnitude of upstream velocity profiles. 

Moreover, the model permitted identification of the observed upstream influence as the 
product of a balance between inertial and gravitational forces, thus distinguishing it from 
the recirculating blocking column produced when viscous and gravitational forces 
predominate. 
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The numerical model was unsuccessful, however, in predicting the propagation of an inter- 
nal solitary wave ahead of the upstream influence. Evidently, to successfully predict the 
generation of such a wave, more work is needed on accurate numerical representation of 
the impulsive start of the plate. 

The model also failed to describe the existence of the reverse lee jet. As mentioned above, 
this could be accomplished by making the plate infinitesimally thin, but the procedure 
introduced singularities which eventually destroyed the numerical solution both upstream 
and downstream. 

CONCLUSIONS 
Summary of Findings 

The significant conclusions of this study are as follows: 

(a) The increase during a series of experimental runs of the gradient-layer thickness 
a' 1 was very nearly a linear function of the number of transits made through 
the layer by the vertical flat plate (see figure 5). 

(b) The observed upstream influence propagated as an internal wave of zero 
frequency whose wave speed depended on the Froude number F = U/(g/a) 1/j . 

(c) The behavior of the observed upstream influence can be adequately predicted by 
a numerical solution of the equations formulated in the section, Theoretical 
Study. This includes a description of the propagation of the upstream influ- 
ence to steadily increasing distances from the obstacle, and the prediction of 
the shape and magnitude of the upstream velocity profiles. 

(d) The Boussinesq approximation can properly be used in the formulation of 
equations to describe the observed upstream influence; moreover, neither mass 
diffusion nor viscosity are important in the upstream flow. 

(e) The upstream influence observed in the present experiments has been identified 
as the result of a balance between inertial and gravitational forces. Furthermore, 
no recirculating flow was observed upstream. The observed flow is thus of a 
different regime from the recirculating blocking column produced in very slow 
flows where viscous and gravitational forces predominate. 

(f) The impulsive start of a vertical flat plate in a region of suddenly varying 
density will produce an internal solitary wave (Davis and Acrivos 1967; 

Benjamin 1967). The numerical flow model used in this study, however, 
fails to predict the generation of this wave. 

(g) A reverse lee jet was observed downstream of the moving plate. The mechanism 
for the formation of the jet is thought to be a buoyancy effect. The present 
numerical flow model with its plate of finite thickness fails to predict the 
existence of the lee jet. 
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Suggestions for Future Research 

The dilemma confronted in attempting to predict numerically the existence of the reverse 
lee jet, without destroying the solution with singularities, suggests that further study is 
needed in the area of numerical flow modeling in the vicinity of mathematically sharp 
edges. 

Another area needing study is the internal structure of the upstream flow when Re = 0(1 ) 
so that the inertial and viscous forces are of the same order. Of interest would be the 
nature of the transition between a recirculating viscous blocking column and a laminar, 
inertial upstream influence. 

Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland May 1973 
502-21-28-01-51 
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NOMENCLATURE 


A overrelaxation constant 

a, c constants used in the mapping from stretched to real coordinates 
a solitary-wave amplitude 

b half-height of the vertical flat plate 

c g internal-wave group velocity 

c° w propagation speed of an infinitesimal, internal solitary wave 

c sw propagation speed of an internal solitary wave of finite amplitude 

c u . propagation speed of the upstream influence 

D mass diffusion coefficient 

F U/(g/a)* , Froude number 

Fj U/(g A p/p H) 1/4 , or ( 1 /Ri) Vl , internal Froude number 

g acceleration of gravity 

H symbol used by Long for total depth of a stratified fluid 
i x-grid index in the finite-difference solution 

j z-grid index in the finite-difference solution 

lc wave number vector of an internal wave 
L b/(Ui>/j3g) 1/3 , dimensionless parameter used by Janowitz 
N (-g/p dp/dz) H , Brunt-Vaisala frequency 

p pressure 

Re U/av, Reynolds number 

Ri (N/Ua) 2 , Richardson number 

s number of iterations through overrelaxation routine 
Sc u/D, Schmidt number 

t time 

t Uat, dimensionless time 

u horizontal fluid velocity with respect to the vertical flat plate 
u u/U, dimensionless horizontal velocity 
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NOMENCLATURE (Continued) 


u dimensionless horizontal velocity with respect to stretched coordinates 
U velocity of the vertical flat plate with respect to the experimenter 

w vertical fluid velocity with respect to the vertical flat plate 

w w/U, dimensionless vertical velocity 

w dimensionless vertical velocity with respect to stretched coordinates 

x, z Cartesian coordinates 

x, jT ax, az, dimensionless Cartesian coordinates 

x, z stretched coordinates 

a a parameter of the density profile; a" 1 is a measure of the gradient-layer thickness 

0 -l/p~ (dp/dz), symbol for constant-density gradient used by those investigators who 

have considered linear density profiles 

At finite-difference time increment 

Ax finite-difference x-grid mesh (stretched coordinates) 

Az finite-difference z-grid mesh (stretched coordinates) 

f z component of vorticity 

17 y component of vorticity 

77 T7/Ua, dimensionless vorticity 

X (g/2ac^ w ) Cn (p c /p u ), dimensionless parameter used by Davis in his study of 
internal solitary waves 

p viscosity 

v p/p, kinematic viscosity 

£ x component of vorticity 

p fluid density 

p fluid density at the horizontal centerline, z = 0 
p p/p, dimensionless density 

p 5 density of the homogeneous salt-water layer 
p u density of the homogeneous fresh-water layer 

a frequency of an internal wave 
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NOMENCLATURE (Continued) 

\p stream function defined by u = d\p/dz, w = -di ///9x 
\p a./ U, dimensionless stream function 

c3 a parameter of the density profile (p 2 - p u )/(p c + p u ) 
uj vorticity (VX^= (£, 77 , f)) 
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APPENDIX 1 


TRANSFORMED GOVERNING EQUATIONS 


The coordinate transformation (19) applied to the governing equations (20) produces the 
following transformed system. The continuity equation becomes 


— + a (1 - x 2 ) • c (1 - z 2 ) 

at 


'a a _ ■ 

— (up) + — (wp) 
3x dz 
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The transformed momentum equation is 
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The velocities are 


d\p _ a^ 

u = — , w = - — 
dz 3x 


(A3) 


and the vorticity expressed as the Laplacian of the stream function is 
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The transformed boundary conditions for t > 0 are 
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The transformed initial conditions are 
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APPENDIX 2 


FINITE-DIFFERENCE FORM OF TRANSFORMED GOVERNING EQUATIONS 


The stream function i p is determined by overrelaxation of (A4) subject to the boundary 
conditions (A5) and (A7). The overrelaxation scheme is given by 


*g +I) = + A ^ij 



(Bl) 


where i and j are the x and z grid indices, respectively, so that is the value computed 

for the stream function at the grid point (i, j) on the (s+ljriz iteration. The value of 1 //. ^ 
is obtained from the finite-difference form of (A4). If the overrelaxation converges, then 
as s increases ' 
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A is an overrelaxation constant whose value affects the speed of convergence but not the 
final result for For this study, A was set .equal to 1.75. The overrelaxation equa- 

tion for the stream function is therefore given by 



= (1 - A) 




-(Ax) 2 ^ 


+ a 2 (l-x . 2 ) 2 + 


- a 2 x. (1 - x. 2 ) Ax O/'j+u - f UJ ) 


(B2) 


+ 


(1-zf) 





- C 


2 v 


Z. (1 - Z ?) 
J v J 7 


(Ax) 2 

(Az) 


(*U» - * 




The boundary conditions on \p at the solid boundary of the obstacle are derived as 
follows. The i// = 0 streamline is coincident with the horizontal centerline of the flow 
field except at the obstacle where it bifurcates and follows the contour of the solid 
boundary. The forward-difference expression for the derivative of a function f at a 
solid boundary is 


~ [- 3f , + - f M ] + 0 [(**>*] 

and the backward-difference expression is 
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(Conte 1965, p. 1 10). The nonslip condition, together with the additional constraint that 
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\p = 0 along the solid boundary, combines with (B3a) and (B3b) to give the following 
boundary conditions for \p. 
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where the obstacle is'represented as a rectangle contained within the grid lines i = M3 to 
i = M4 and j = N3 to j = N4. The u and w velocity fields, as given by (A3), become 
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The continuity equation (Al) can be written in finite-difference form as 
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where the unprimed variables represent values at t = t and p! . is the new value of the 
density at grid point (i, j) and at t = t + At. As mentioned in the section, Theoretical 
Study, the nonlinear terms [3/33T(up)]. . and [d/dz~(c op)]. . were evaluated by the 
“special three-point noncentral difference method” of Torrance and Rockett (1969). 


A derivation analogous to that for the boundary conditions on \]j gives the following 
boundary conditions for p. 
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The momentum equation, which is used to evaluate the change of free stream vorticity 
with time, can be written in finite-difference form as 
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Finally, the vorticity at the solid boundaries is evaluated from i? = V 2 \jj by taking note of 
the boundary conditions on \ p and using the forward- and backward-difference expressions 
for a second derivative at a boundary. These expressions are (Conte 1965, p. 110) 
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APPENDIX 3 


COMPUTER PROGRAM 

The computer program used to make the theoretical flow predictions for this study is 
based on the finite difference equations of Appendix 2. It is written in FORTRAN-IV 
for use on an IBM 360 computer. It uses double precision, and for the present 81 by 41 
grid mesh, requires 620 K bytes of storage and 30 to 40 minutes of central processing 
unit (CPU) time to model a complete experimental run. Most of the CPU time is spent in 
the overrelaxation step of the computation. The core and CPU time requirements are 
considerably reduced for coarser grid meshes, but not without some sacrifice in accuracy. 

INPUT-DATA FORMAT 

Cards 1, 2, 3, and 4 are title cards. A descriptive title can be punched on these cards in 
columns 1 to 80. 

Card 5 contains 12 floating-point data inputs in fields of six columns each: 

Col. 1 - 6 DX = Ax = finite-difference x-grid mesh (“stretched” coordinates) 

7-12 DZ = Az = finite-difference z-grid mesh (“stretched” coordinates) 

13-18 DT = At = dimensionless finite difference time increment 

19-24 T = time at start of computation 

25 - 30 F = Froude number = U/(g/a) yz 

31-36 Re = Reynolds number = U/w 

37 - 42 EPS = (Schmidt number)' 1 = D/v 

43 - 48 DAVG = p = fluid density at z = 0 

49 - 54 OMEG = cc = a parameter of the density profile 

55 - 60 ALPHA = a = a parameter of the density profile 

61-66 UBODY = U = velocity of the vertical flat plate with respect to the 
experimenter 

67-72 A = overrelaxation constant 

Card 6 contains nine fixed point data inputs, right justified, in fields of six columns each: 

Col. 1-6 LL = total number of time incrementations to be made during the 
computation (should be a multiple of 40) 
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7-12 MX = total number of x-grid lines; must be less than or equal to the 
number of storage locations assigned to X(I) in the DIMENSION 
statement 

13 - 18 ' M2NEG = downstream thickness of the obstacle in terms of the number 
of downstream (x < 0) grid lines it contains 

19-24 M2POS = upstream thickness of the obstacle in terms of the number of 
upstream (x > 0) grid lines it contains 

25- 30 NZ = total number of z grid lines; must be less than or equal to the 
number of storage locations assigned to Z(J) in the DIMENSION 
statement 

31-36 MLT1 = 0, if Boussinesq approximation is to be made 
= 1 , otherwise 

37 - 42 MLT2 = 0, if Boussinesq approximation is to be made 
= 1 , otherwise 

43 - 48 MLT3 = 0, if flow is to be assumed inviscid 
= 1 , otherwise 

49 - 54 MLT4 = 0, if flow is to be assumed nondiffusive 
= 1 , otherwise 

Card 7 contains 12 fixed point data inputs, right justified, in fields of six columns each. 
These input parameters indicate those results which the user has selected to be printed out. 
If a parameter is zero, the results corresponding to this parameter are not printed out. If 
the printout is desired, a nonzero integer must be assigned to the parameter. The following 
list gives each parameter in the order in which it is read in and the result prints out to which 
it corresponds: 

Col. 1 - 6 IPOK - number of iterations K through SUBROUTINE RELAX 
7-12 IPOP - stream function i p both initial and at later time 
13-18 IPOPC - stream function contour elevations 
19-24 IPODI - initial density profile at a typical cross-section 
25 - 30 IPODC - density contour elevations 
31-36 IPODN - nondimensional density at later time 
37-42 IPOD - dimensional density at later time 
43 - 48 IPOV - vorticity field 

49 - 54 IPOUO - nondimensional u-velocity with respect to the obstacle 
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55 - 60 IPOWO - nondimensional w-velocity with respect to the obstacle 
61-66 IPOUE - dimensional u-velocity with respect to the experimenter 

67- 72 IPOWE - dimensional w-velocity with respect to the experimenter 

COMPUTER PROGRAM LISTING 

The following pages present the program listing. 
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WRITE (6*1090) A*LL*MX*M2NEG*M2P0S»NZ MAN 32 

WRITE (6*1100) MLT1*MLT2*MLT3*MLT4 MAN 33 

IF (MLT1) 10*20*30 MAN 34 

10 WRITE (6*1110) MAN 35 

GO TO 1010 MAN 36 



20 WRITE (6*1120) MAN 37 

60 TO 40 ' MAN 38 

30 WRITE (6*1130) MAN 39 

40 IF (MLT2) 50*60*70 MAN 40 
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DO 170 I = 1 » MX MAN 104 
X ( I ) =DX*UFL 0 AT ( I-M 2 ) MAN 105 
XSEC 2 (I>=D 1 -X<I )»*2 MAN 106 
Y=X(I) MAN 107 
XR(I)=(D 1 /CA)^ARCTNH(Y) MAN 108 



XR(1)=(-1.0DS)*ALPHA MAN 109 

XR( MX )=(1#0D6)* ALPHA MAN 110 

XROIM=OIMENSIONAL FORM OF XR MAN 111 

XRDIM(I)=XR(I) /ALPHA MAN 112 
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2J) ) MAN 380 

CALL NONLIN (DUVDX »UAV » VT 1 * DX » DWVDZ » WAV *VT 1 *DZ * I* J) MAN 381 

CALL DELSQR < DSQV , CA»X » VT 1 ,DX» ALF» Z » DZ » I * J) MAN 382 
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IF (IPOV.EQ.O) GO TO 780 MAN 447 
WRITE (6* 1590 ) MAN 448 
WRITE (6*1600) T »UIMT MAN 449 
WRITE (6*1640) MAN 450 
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